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Abstract 

The Spectral Volume (SV) method is extended to the 2D Euler equations. The focus of this paper 
• a thp nprformance of the SV method on multidimensional non-linear systems. 

Wementtion details including total vanation diminishing (TVD) and total variation bounded 
(TVB) limiters are presented. Solutions with both smooth features and discontinuities are uti ize 
to demonstrate the overall capability of the SV method. 

1. Introduction 

We continue the development of the Spectral (Finite) Volume (SV) method for hyperbolic 
Z : s ~)lws unstructured gods. Ultimately, the SV 

volume method which has been under development for several decades, and has become tne 
stateof-the-art for the numerical solution of hyperbolic conservation laws. Similar to he 
Godunov method, the SV method has two key components One is data “““truaiom and 

volumes are partitioned in a geometrically similar manner in a simplex, and thus a sing 
reconstruction is obtained. Next, the DOFs are updated to high^rder accuracy^singft^usud 
Godunov method. Numerical tests with scalar conservation laws in both ID and 2D an 
Systems hTve verified that the SV method is indeed highly accurate, conservative, and 

geometrically flexible. 

2. Spectral Volume Method for the 2D Euler Equations 

The unsteady 2D Euler equation in conservative form can be wntten as 
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where Q is the vector of conserved variables, E and F are the inviscid flux vectors. Assume that 
we solve (2.1) in the computational domain G subject to proper initial and boundaiy conditions. 
The domain Q is discretized into N non-overlapping triangular spectral volumes (S s) 

Q = U 5,. 

Given a desired order of accuracy k for (2.1), each spectral volume S, is then partitioned into m = 
k(k+l)/2 control volumes (CVs), and th ej-th CV of S, is then denoted by C tJ . Let q denote any 
of the conservative variables in Q. The cell-averaged conservative variable q at time t for control 
volume Cij is defined as 


j q(x, y,t)dxdy 
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where V Ui is the volume of C„ ; . Given the cell-averaged conservative variables for all the CVs 
in Si, a polynomial p t (x,y)e P k 1 (the space of polynomials of degree at most k - 1 ) can be 
reconstmcted such that it is a k-th order accurate approximation to the function q(x,y) inside So 

P i (x,y) = q(x,y) + 0(h k ), (x,y)eS t , i = l,...,N. (•) 

This reconstruction can be solved analytically by satisfying the following conditions: 
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This polynomial Pi (x,y) is the k-th order approximation we are looking for as long as the function 
smooth in L region covered by S„ The reconstmciion can be more convemently 

expressed as 

P i (x,y) = ^L j (x,y)q iJ , ^ 

where Lj(x,y)e P k ~ l are the "shape" functions which satisfy 

J L m (x,y)dxdy 
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The high-order reconstruction is then used to generate high-order updates for the cell-averaged 
state variable on the CVs. Integrating (2.1) in C u , we obtain the following integral equation for 

the CV-averaged mean _ 

%- + T7-L|(/*" )<M = 0 ' 
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: .j '■=' A, 
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where O is the vector of the CV-averaged conservative variables in C u ,f- (E, F), K is the 
number of faces in Cy, and A, represents the r-th face of Cy. The surface integration on each face 
can be performed with a k-th order accurate Gauss quadrature formula, i.e. 

| (/ • n)dA = ^w rn f(Q(x rq y rq )) •n r A r + 0(A r h k ), 


( 2 . 8 ) 


where J - integerKk+mi is the number of quadrature points on the r-Ih face and wv, are the 
Gauss quadrature weights/ (x rq , y r ,> are the Gauss quadrature points h is the ■ ">™ edge 
length of all the CVs. Time t is omitted whenever there is no confusion. / 
following identity exists: 

(2.9) 


£}(/ 7 *n)dA = 0. 
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Therefore, we will gain an extra order of accuracy if we sum up the surface integrals for the 
faces of Cjj, i.e., 

V f (/ • n)dA = XXvv r9 /(G(x r<7i y r( ,)) • «A + 0{Afi k * x ). 

r=l 1 r= 1 -J=l 

Since 0(V,) = 0(A r h) , we therefore have 

(/ ,n)dA = 2- 1 1 + >• 

With the polynomial distribution on each SV, the state variable is most likely discontinuous 
across the^ SV boundaries. Therefore, the flux integration involves two discontinuous state 
varies iust to the left and nght of a face of the SV boundary. This flux integration ,s earned 
out using an exact Riemann solver or one of the Lipschitz continuous approximate Riemann 
solvers or flux splitting procedures, i.e., 

/(Q( X r^y r q)) 0n r ~fmm(QL( X rq>yr q )’QA X rg^yr q )’ n r^ 

where Q h and Q R are the vector of conservative variables just to the left and nght of a face. 
Substituting (2.12) into (2.11), we obtain 

= (2-13) 

It is the Riemann solver which introduces the "upwinding", and dissipation into the SV method 
such ftat the SV method is no. only h.gh-order accurate, but also stable. In this paper, we employ 
Roe* and* local Lax-Friedrichs flux formula. Finally we obtain the following sem,-d,screte SV 

scheme _ 

<2 ' 14) 

For time integration, we use the 'third-order TVD Runge-Kutta scheme from. 

3, Multi-Dimensional TVD and TVB Limiters 

For the non-linear Euler equations, it is necessary to perform data limiting to maintain stabi lty i 
£ soMon contain, discontinuities. There are two different ways 

system setting. One way is to apply a limiter to each characteristic vana • P? 

a limiter to each component of the vector of the conservative variables. The former 
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property of naturally degenerating to the scalar case if the hyperbolic system is linear, but 
latter is much more efficient. In this paper, we choose the component-wise approach because 
its efficiency. To this end, we first establish the following numerical monotonicity criterion tor 

each control volume (3 1) 


q- < q!j (x rq ,y yq )<q- 


where q , ™ n and q ™ “ are the minimum and maximum cell-averaged solutions among all its 


(3.2) 


neighboring CVs sharing a face with Cq , i.e., 

= max(^ r ) 

q!j = rmn(q tJ ,mmq ljr ), 
and a (x y ) is the reconstructed solution at any of the quadrature points. If (3.1) is strictly 

enforced, the resultant numerical scheme for the scalar case is TVD. However, it is well known 
that TVD schemes are locally first-order at extrema, and may pollute the globa accuracy o 
solution. In order to maintain the order of accuracy away from discontinuities the TVB idea 
followed here, i.e., small oscillations are allowed in the solution. If we express the reconstruction 
for the quadrature points in the following form 

M rq = PX X r^yr q )-^,P 

then no data limiting is necessary if 


H, S4MX- 


where A =U,-r,| is the distance from the CV centroid to the quadrature point. In other 

words, no data limiting is necessary if (3.3) is satisfied, even if (3.1) is not. Usually the constant 
M, is chosen to be the maximum second derivative of the solution. However M, is a user chosen 
parameter if there is a discontinuity in the solution. Note that a different U, should be used for a 
different conservative variable. In this paper, M, is scaled according to the minimum and 
maximum of the component, i.e., 

M q =M(g nlax ; 

where M is a constant independent of the component, and q max and q min are the maximum and 
minimum of the solution q over the computational domain. If (3.3) is violated for any quadrature 
point, then it is assumed that the CV is close to a discontinuity, and the solution in the CV 

locally linear, i.e., ~ 

= ;+V«u*fr- r u ), VreC,.,. (3.5) 

The magnitude of the solution gradient is maximized subject to the monotonicity condition given 
in (3.1). The original high-order reconstruction in the CV is used to compute a initial guess of the 

gradient, i.e.. 




d P< 


\ 


dx ’ dy 

\ ' -i ,j 

This gradient may not satisfy (3.1). Therefore it is limited by multiplying a scalar <pe [0, 1] so 

that the following solution satisfy (3. 1) 

q IJ (x,y) = q i ' J +<pVq u '( r - r u)- ^ 

In the case of M = 0, the limiter becomes TVD. The availability of cell-averaged data on the CVs 
inside a SV makes this CV-based data limiting possible, whereas in the DG method, one can 
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only do an element based data limiting. Due to the increased local resolution, the SV method was 
shown to have better resolutions for discontinuities than the DG method. This advantage s ou 
also carry over to the multi-dimensional systems case. 

4. Numerical Tests 

* unit wide, with a step of 0, units ^ ,ca t e d at 0 6 
units from the tunnel inlet. The initial condition is a Mach 3 right-going uniform flow. In viscid 
wall boulary conditions (reflective) are used for tunnel wall boundaries, while inflow and 
outflow boundary conditions are used at the inlet and exit of wind tunnel It is well taiowntha^ 
the comer of the step is a singularity, and can lead to a spurious Mach stem at the downstr 
bottonTwall, and anerroneous entropy layer a. the bottom wall. In the presen study no^ a 
treatments were used for the singularity to see how the singularity affects 
solutions. 

Four different gnds with various gnd resolution near the singular comer point 1« 

The first grid is nearly uniform with a mesh size of h = 0.025. The other three gn 
Successively refined by halving the grid size near the comer. The four gnds are shown in Figure 
“^computed density distributions at t = 4 with the second order SV scheme and a TVB 
limiter (M = 10) is shown in Figure 2, and the corresponding third-order results are displayed t 
Figure 3 First it is evident that the spurious Mach stem weakened considerably with increased 
Ji resolution near the comer point in both the second-order 

The entropy layer downstream the shock reflection point also weakened. Note that the third order 
results have higher resolutions for the shock waves, as expected. 

this problem is chosen to be [0, 4] x [0, 1], The reflecting wall lies 
I the boSom of the computational domain starting from *=1/6. Initially a nght-movtng ; Mach ,10 
shock is positioned at *=1/6, y=0 and makes a 60" angle with the For ^the bottom 

boundary the exact post-shock condition is imposed for the region from x-0 to x-1/6 and a solid 
wall boundary condition is used for the rest. For the top boundary of the computational domain 
the solution is set to describe the exact motion of the Mach 10 shock. The left boundary 
the exact post-shock condition, while the right boundary is set as outflow boundary. The results 
“Tshown. Figure 4 shows .he density contours computed wHb . Ito second ^SV scheme 
and a TVB limiter (M = 10) on the coarse (Ax=Ay=l/60), medium (Ax=Ay=l/120) and fine 
(Ax=Ay=l/240) grids. The corresponding results for the third SV scheme are shown in 'gmc - 
if* obvious that the third order SV scheme captures the shock waves and the complex flow 
structure under the triple Mach stem better than the second-order SV scheme. 

5. Conclusions 

The Spectral Volume method has been successfully extended to two-dimensional hyperbolic 
lystems of conservation laws. The extension to the system has been perfomted in a component- 
wise manner Both TVD and TVB limiters are also implemented in a component-wise manne 
This approach is very efficient, and has been shown to produce good numerical results 

smooth and non-smooth flow problems. 
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Figure 1. Four computational grids with different degrees of refinement near the singular comer 

point 
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Figure 4 Density contours computed with the second order SV schemes and a TVB limiter 
6 (M = 10) 


9 




1 




